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Abstract. The class of norm-dependent Random Matrix Ensembles is studied in 
the presence of an external field. The probability density in those ensembles depends 
on the trace of the squared random matrices, but is otherwise arbitrary. An exact 
mapping to superspace is performed. A transformation formula is derived which gives 
the probability density in superspace as a single integral over the probability density 
in ordinary space. This is done for orthogonal, unitary and symplectic symmetry. In 
the case of unitary symmetry, some explicit results for the correlation functions are 
derived. 
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1. Introduction 

Supersymmetry is a prominent and widely used tool in studying disordered systems 
and systems that can be modeled by random matrices, see Refs. jll El El IH IHl E] • The 
method was developed for Gaussian probability densities, a review and a discussion 
of the mathematical justification was recently given in Ref. [7j. This restriction to 
Gaussian probability densities is no shortcoming if one is exclusively interested in 
calculating correlations on the local scale of the mean level spacing. This is due to 
local universality [HI E] ■ Probability densities which do not introduce scales competing 
with the mean level spacing yield correlations which are one the local scale identical to 
the ones resulting from Gaussian probability densities, see a review in Ref. When 
studying matrix models in high-energy physics one is not interested in the local scale. 
Another universality in the leading asymptotics of the matrix dimension was found ^U] 
for the correlation functions on large scales. 

Nevertheless, restriction to Gaussian probability densities does not always suffice. 
First, the one-point functions obviously depend on the specific form of the probability 
densities, because they are not measured on local scales. Such level densities are 
important, for example for applications in high energy physics ^T], but also in 
finance ^2)- Second, to distinguish certain directions in matrix space, one adds 
an external field to the random matrix, and one often averages over the matrices 
representing the external field. The local correlations now change and depend sensitively 
on the root-mean square matrix element of the external field divided by the local mean 
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level spacing. Examples are the crossover transitions, see the review in Refs. jH E] • In 
such a situation, scales competing with the local mean level spacing might occur which 
can lead to a deviation from universal features, such that the crossover transitions would 
differ for different probability densities. Third, non-Gaussian probability densities and 
their non-universal features on special scales have always been of interest in conceptual 
studies and for considerations in general statistical mechanics, we mention the bound- 
trace and the fixed-trace ensembles JH] and the recently introduced ensembles deriving 
from a non-extensive entropy principle . 

In the present contribution, we show that the supersymmetry method can be 
extended to random matrix models with non-Gaussian probability densities. In the 
context of universality, asymptotic considerations for infinite level number have already 
been combined with super symmetric techniques for non-Gaussian probability densities 
in Ref. ^Hl- Here, however, we aim at an exact discussion. We focus on the large 
class of norm-dependent random matrix ensembles which depend through an arbitrary 
functional form on the trace of the squared random matrices. Recently, a general 
construction of these ensembles was given in Ref. J7] . For the reasons just mentioned, we 
include an external field. We have two goals. First, we want to deliver the conceptually 
important insight that supersymmetry is by no means restricted to Gaussian probability 
densities. Second, we want to provide a series of explicit and practically relevant 
formulae for the correlation functions. Here, the application of supersymmetry yields 
particularly handy results in the presence of an external field. 

The question whether or not norm-dependent ensembles can be formulated exactly 
in a supersymmetric framework was also discussed by F. Kalisch. Although his approach 
was quite different from the one to be presented here, it would have been likely to 
produce equivalent results for the case without external field. Unfortunately, F. Kalisch 
left academia and his findings are unpublished. 

The paper is organized as follows. We formulate the problem in Section El thereby 
also introducing our notations and conventions. The supersymmetric representation of 
the norm-dependent ensembles is constructed in Section 01 In Section 01 we discuss a 
series of examples. Explicit results for the correlation functions are given in Section 
We summarize and conclude in Section (HI 

2. Formulation of the Problem 

In Section \'2.1\ we set up the generating function in the presence of an external field. 
As we need to refer to the Gaussian case, we briefly sketch it in Section IT^ We discuss 
norm-dependent ensembles and pose the problem in Section 12.31 

2.1. Generating Function in the Presence of an External Field 

The three symmetry classes of N x N random matrices H are labeled by the Dyson 
index jS. In the orthogonal class, H is real symmetric {f5 = 1) and in the unitary class, H 
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is Hermitean (/3 = 2). In the symplectic class, H is self-dual {(3 = 4) and the entries of 
H are 2x2 quaternions. The eigenvalues of H are doubly degenerate in the symplectic 
class. We notice that such a matrix has 

^ = N + ^NiN-l) (1) 

independent matrix elements. The quantity fi is often referred to as the number of 
degrees of freedom. A normalized probability density function P^^\H) defines — 
together with the symmetry class — the random matrix ensemble. We add a fixed 
external field represented by a matrix Hq which, without loss of generality, can be 
assumed to be diagonal. In the symplectic class, it has dimension 2N x 2N . Thus, 
we are interested in a system described by Hq + aH where a measures the relative 
strength. The k level correlation function R^if\xi, . . . ,Xk,a, Hq) is the probability 
density to find k eigenvalues of Ho+aH at positions xi, . . . ,Xk- The correlation functions 
i?[f^(xi, . . . , Xfc, a, Hq) are technically easier to handle. They include real and imaginary 
parts of the propagator, while the R^f^\xi, . . . ,Xk,a, Hq) are only the correlations of 
the imaginary parts. The latter can easily be constructed from the former. We use the 
conventions of Refs. ^11021 EO]- For arbitrary P^^\H), the correlation function can be 
written as the derivative 



^ 1 0^ 



(2) 



of a generating function 



zfix d[H] p<^\H) n i TuHTH'^'''ln )'' 

J pJi \ det{Ho + aH -Xp + Jp) J 

with respect to source variables Jp, p = 1, . . . , fc. Here, we define 7 = lif/9 = l,2 and 
7 = 2 if /3 = 4, moreover we introduce the diagonal matrices x = diag {xi,xi, . . . ,Xk, Xk) 
and X = diag (— Ji, + Ji, . . . , —Jk, +Jk)- The volume element d[H] is simply the product 
of the differentials of all independent matrix elements. For complex variables, we use 
the differentials of real and imaginary part. 

2.2. Gaussian Random Matrix Ensembles 

In the Gaussian case, the normalized probability density function with variance 2f^//3 
reads 

P^'^-iH)^^,(^)f'\.,(-l,^H^) (4) 

with the number ^ of degrees of freedom given in Eq. (^. To properly account for the 
degeneracies in the symplectic class, we define 

if /3 = 1,2 

Tr = <( 1 . (5) 

^ tr if /5 = 4 ^ ' 
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(6) 



The probability density (0]) and the symmetry class define the Gaussian Orthogonal, 
Unitary and Symplectic Ensemble GOE, GUE and GSE for /5 = 1,2,4, respectively. 
The generating function 

Zf''\x + J, 2vy(3) = J d[H] P(^^)(i7) 

'det{Ho + aH - 
det{Ho + aH - 

of this case has an exact representation as integral in superspace, 

+ J, 2vyf3) = J d[a] Q^^''\(t) 

detg -^/'^ {{aa - x - J) ® I^n + l(k ® Hq) . (7) 

The crucial feature of supersymmetry is the drastic reduction in the number of degrees 
of freedom. This is borne out in the dimension of the matrix a. It is a 2/c x 2k Hermitean 
supermatrix for (3 = 2 and a. 4k x 4k Hermitean supermatrix with additional symmetries 
for /3 = 1, 4 1^, we use the conventions of Refs. pT| l22j. The parameter ( = 2 for (5 = 2 
and C = 4 for /? = 1, 4 is defined accordingly. We write 1m for the M x M unit matrix. 
Thus, the expression ((7j) contains the unit matrices l^v and l^k for /5 = 1, Iat and 
for /5 = 2 and I2N and l^k for /? = 4. Again, the volume element d[a] is the product of 
the differentials of all independent variables. For the complex anticommuting variables, 
we use the differentials of the variable and of its complex conjugate. The probability 
density in superspace 



QiGp) = ^W) ( -£^trg a' ) (8) 



is a normalized Gaussian as well. Importantly, the normalization constants 

\ 2'=(4/c-3)/2 if (3=1,4 

depend only on the dimension. In contrast to the ordinary case (@)), they doe not contain 
the variance 2f ^/ (3. 

The result ((Tj) has a remarkable property. The superdeterminant comprises a sum 
of two terms which are both direct products. The first term is a direct product of 
supermatrices with a unit matrix and ordinary space, and vice versa in the second term. 
Most conveniently, this decouples to some extent the random matrix ensemble, i.e. the 
matrix a from the external field Hq. This feature, which is typical for the supersymmetry 
method, was already very helpful for an exact calculation of the transition from Poisson 
regularity to the GUE in Refs. 1^ . Furthermore, it also made possible some 
asymptotic evaluation [23 121] of the correlations on the local scale for large coupling 
a/D where D is the mean level spacing. 

2.3. Posing the Problem for Norm-dependent Ensembles 



In analogy to the scalar product for vectors, one introduces a scalar product Tr HK for 
two matrices H and K with the same symmetries. This is then used to define the norm 
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of a matrix by 



II H 11= VTtH'^ , (10) 

corresponding to the length of a vector. The class of norm-dependent ensembles has a 
probability density of the form 

P(^)(i7) = p(^^)(Tri72) , (11) 

where P^^^\u) is function of the norm || H \\ or, equivalently of u = TtH^. Of course, 
piTf3) ]-^g^g chosen such that P*-^-* (i/ ) is positive semi-definite and fulfills all the 
necessary convergence requirements. According to the symmetries, there are Norm- 
dependent Orthogonal, Unitary and Symplectic Ensembles for (3 = 1,2,4. We denote 
them TOE, TUE and TSE, respectively. We show in Appendix Aj that the u-th moment 
of the probability density can be expressed, if it exists, as the single integral 

This includes the normalization by setting Mq = 1 for z/ = 0. 

Many ensembles fall into the norm-dependent class. Obviously, the Gaussian 
Ensembles are found by setting P^'^^\u) ~ exp(— /5M/4f ^), which is an exponential 
function, not a Gaussian. Non-trivial examples are the fixed-trace and the bound-trace 
ensembles ^S]- An important subclass of norm-dependent ensembles is derived from a 
non-extensive entropy principle . It comprises a variety of interesting cases which 

are found by considering limits of certain parameter. A rather general construction of 
norm-dependent ensembles using a single-valued spread function is given in Ref. |17j . 
We return to this point. 

We ask the following questions. Can we express the generating function for the 
norm-dependent ensembles TOE, TUE and TSE given by 

Zp)(a; + J)= I d[if]p™ (Tr/J2) 

' det{Ho + aH -Xp- Jpj V ^^^-^ 
det(ifo + ctH — Xp + Jp) J 
as integral in superspace? — Can we construct the supersymmetric analog of the 
probability density P^'^'^\Tt H^) ? — The answers are in the affirmative. We will 
derive the exact representation 

ZP^x + J)= j d[a] Q(^^)(trga2) 

detg-^/2^((aa-x- J)® V + lcfc®^o) , (14) 
where the supermatrices a are defined as above and where the probability density 
Q(-^^)(trgcr^) is also norm-dependent, but now in superspace. 

Importantly, the direct product structure implying the decoupling of the random 
matrix ensemble from the external field Hq is also present here for all TOE, TUE and 
TSE. This extends the discussion in Section 12.21 for the Gaussian ensembles. 
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3. Supersymmetric Representation 

In Sectioning we derive the supersymmetric representation by using Fourier integrals. 
We present the transformation formulae for the probability densities in Section 13.21 
The connection to the spread function is discussed in Section 13. 3| which also contains 
an alternative derivation of the transformation formulae. 



3.1. Derivation Using Fourier Integrals 

The norm || H \\ is non-negative and we have u = Ti > 0. Thus, P^'^^\u) is only 
defined on the positive u axis. When introducing the Fourier integral over the entire 
axis, we have to set P^^^\u) = for u < 0, such that 

CO 

{iyu)du (15) 

v27r J 
^ 

is the Fourier transform with the inversion 

+ CO 

piTP)^u) = [ p^^^\y) exp {-lyu) dy . (16) 

v27r J 

— CO 

We add a small imaginary increment to the Fourier variable, y~ = y — ie and insert 
insert Eq. ()16p with u = Trif^, i.e. the integral 

+ 00 

pim H^) = —= f (y) exp (-ly- Tr H^) dy . (17) 



27r 

-oo 

into the generating function (fT^ . We thereby rediscover the Gaussian case (0) with 
the variance l/i2y~. The integrals over H can now be done as Gaussian integrals, 
the complex variance l/i2y^ does not cause a problem. Even without the imaginary 
increment, they exist as Fresnel integrals. The imaginary increment makes standard 
Gaussian integrals out of them, but this is not the motivation for it. We need 
the imaginary increment later on. The important difference to the Gaussian case of 
Section 1221 is the fact that the Gaussian exp {—iy~ Tr H^) comes without normalization 
constant for the H integration. Hence, when inserting Eq. ()17|) into Eq. (jH)) we obtain the 
inverse of the normalization constant as an y dependent factor in the Fourier integral. 



ZP)(X + J) = ^ / dyp^^^\y) 
v27r 



2^/^(^^j Zr\x + J,l/^2y-) (18) 

with Zl^^\x + J,l/z2|/-) given in Eq. (jH)). We now employ the supersymmetric 
representation ((Zj) and find 

Zr)(x + J) = -i= / dyp(^^\y)2^/^l- 



Norm- dependent Random Matrix Ensembles in External Field and Supersymmetry 7 

J d[a] c*-^^ exp (^—iy^trga'^^ 

detg-^/2^((aa-x- J) ® V + Icfc^o) • (19) 

Hence, by interchanging the integrations, we arrive at the desired Eq. ()14|1 . where the 
probabihty density in superspace 

^ c«^.^ i^y" J/^-^ r^g;f:f\ . (20) 

— oo 

is the inverse Fourier integral with an additional power (iy~)^^'^ in the denominator. 
3.2. Transformation Formulae 

We set w = trg cr^ and plug the Fourier integral (fTH|) into Eq. (f^ . 

-oo ^ ^ ' 

The y integral converges because of the imaginary increment and can be done in a 
standard way jSS!- Apart from factors, it yields 6(m — w){u — wY^'^~^. We thus arrive 
at the transformation formula 

= TTT^ (I) h'^'' + w)u^l--Mu . (22) 







This result allows one to calculate, by a single integration, the probability density in 
superspace for any norm-dependent ensemble TOE, TUE and TSE. We notice that the 
Fourier integral (fTS|) has to exist, but, importantly, its explicit knowledge is not needed to 
obtain the probability density Q^'^^^w) in superspace. Interestingly, the transformation 
formula can be inverted. For even number of degrees of freedom /i, iterated integration 
by parts yields 

This inversion is likely to be correct even for odd /i if the theory of fractional derivatives 
is applied. 

From a conceptual viewpoint, the pair of transformation formulae ()22|1 and ()2Hj] 
states the main result of this contribution. The power of supersymmetry lies in the 
drastic reduction of the degrees of freedom. The mechanism of how this happens 
was previously only known in the Gaussian case. The transformation formulae ()22j) 
and ()23p considerably generalize that. A particularly interesting interpretation follows 
from formula ^I^ . The probability densities P^'^^\u) and Q^^^\u) formally coincide 
for = 0, i.e. in zero dimensions, = 0. This is already visible in the Gaussian case. 
Apart from the variance independent normalization c^^^ , the Gaussian (jHj) in superspace 
with u = trg 0"^ indeed follows from the Gaussian Q in ordinary space with m = Tr 
by simply setting = 0. 
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Another interesting observation results from putting w = in the transformation 
formula ()22j) and then using the normalization of P^'^^'>{u) which can be read off from 
Eq. ()12|1 for z/ = 0. One has Q'^'^^\Q) = In other words, the normalization of 

the probability density in ordinary space corresponds to the value of the probability 
density in superspace at = 0. This reflects the Efetov-Wegner-Parisi-Sourlas 
theorem ^1 |23 1211 I2Z]; referred to as Rothstein theorem in mathematics |2H1- It 
implies that the normalization integral for a function such as ours which only depends 
on invariants reads 

l=/g(^^ntrgaVM = ^g(^^)(0). (24) 

This phenomenon exclusively occurs in superspace is due to a subtle mutual cancellation 
of singularities. Hence, it is reassuring to see that the normalization of the probability 
density in ordinary space leads — via the Efetov-Wegner-Parisi-Sourlas theorem — to 
the normalization of the probability density in superspace. 



3.3. Connection to the Spread Function 

A rather general construction of norm-dependent ensembles was given by Muttalib and 
Klauder ^\ for the unitary case. It can be generalized to all three symmetry classes in 
a straightforward manner. The probability density 

P™(Tr/f^) = ]f^^\t) ^ (|^)'''exp (-^Tr/7^) dt (25) 

is expressed as an integral involving a normalized Gaussian with a real variance 2t/l3. 
The quantity /^^^^(t) is referred to as spread function. As seen form Eq. (j^Hl) it is 
normalized, 

oo 

/(^/3)(t)dt = 1 . (26) 

We insert the integral into the generating function (fT!?|) and find in steps analogous 
to the ones in Section 13.11 



oo 

(x + J) = j dt /(^^) it) (x + J, 2t//3) 



oo , 

2 



dtf^^^Kt) j d[a] c^^^exp (-£trg(T^ 
^ ' 

deig-^'''^ {{aa -x~ J)®l^N + hkHQ) ■ (27) 

This yields again Eq. where the probability density in superspace now reads 

Q™(trga2) = J/(^^)(t) c(^) exp (-^trga^^j dt . (28) 
^ ' 

Comparing Eqs. (^3]) and ()28|) one sees that the probability densities are in ordinary and 
in superspace given as integrals over the spread function times a normalized Gaussian. 
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Moreover, we notice that the variable t in the variance 2t/(3 has the meaning of a 
diffusion time. In ordinary space, the diffusion is Dyson's Brownian Motion [211 IHUj . 
It has a fully fledged analog in superspace 120] with the same diffusion time. Thus, 
the TOE, TUE and TSE are, in ordinary and in superspace, ensembles constructed as 
averages involving the diffusion time. 

The transformation formulae (j^^ and (j^H)) are easily re-derived from Eqs. 
and (|^H|). We emphasize that only the existence, but not the precise knowledge of the 
spread function is needed to calculate the probability density in superspace. Those 
readers might appreciate the alternative derivation of the transformation formulae by 
means of the spread function who did not feel comfortable with our treatment of the 
singularities in the Fourier integrals of Section 13.11 

4. Some Specific Examples 

To gain insight into how the transformation formulae work, we calculate the probability 
densities in superspace for a variety of examples. To acquire some first experience, 
we revisit the Gaussian ensembles in Section 14.11 We discuss, for all three symmetry 
classes f3 = 1, 2, 4 the bound trace, the fixed trace, the Gauss-monomial and the Gauss- 
quartic ensembles in Sections 14.21 to 11751 For the probability densities in ordinary space 
of these examples, we introduce constants Oq, Oi and a2 which are always assumed to 
be real and positive. Using Eq. (fT^ . they can be expressed in terms of the moments 
M^l^\ In particular, the overall normalization constant can be fixed with Eq. ()12p 
for 1/ = 0. However, we rather use the relation Q^'^^\Q) = c^^^^ which is according to 
Section equivalent to the normalization of the probability density in ordinary space. 
In Section ^Gl we discuss the ensembles derived from an non-extensive entropy principle. 
We always write u = Tr and w = trgcr^. 

4.I. Revisiting the Gaussian Ensembles 

Inserting the Gaussian into the transformation formula ()22|1 . we find 



can be also be used to infer the spread function, which is a 5 function. The integral ()25p 
has to be interpreted as a proper Cauchy integral. 




(29) 



which gives the Gaussian (jH)). The Fourier transform 




(30) 
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4-2. Bound Trace Ensembles 

The probability density sets a cutoff for tlie norm of tlie random matrices according 

to [in! 

p(s^/5)(^,) = aoe(ai-u) . (31) 
Tlie transformation formula (j^^ yields 

Q^^^^\w) = M^^^^e{a,-u^) , (32) 

Oi 

which is a bound trace ensemble as well, but now in superspace and multiplied with a 
polynomial factor. 

4-3. Fixed Trace Ensembles 

The probability density fixes the norm of the random matrices such that 1^ E21 

P^-^^'P^u) = ao6 {ai - u) . (33) 
With the transformation formula we find 

Q^^^^\w) = M'^'^^^^T' q (ai - w) , (34) 



which is, once more, a bound trace ensemble of the form ()32|) . We notice that the 
exponent in the polynomial factor is /i/2 — 1 compared with /i/2 in Eq. (j32|) . This 
simply reflects that the probability density is the derivative of the probability 
density (jHT|) . We mention that fixed trace ensembles do not seem to exist in superspace, 
at least not in a simple-minded interpretation. This is so, because the normalization 
requirement Q'^'^^\Q) = c^^-' can hardly be fulfilled if Q^'^^\w) includes a S function. 

4.4- Gauss-Monomial Ensembles 

The probability densities in superspace derived in the previous examples tend to have 
remarkable similarity to the ones in ordinary space. This seems to be a fairly robust 
phenomenon. To illustrate it further, we introduce ensembles comprising a Gaussian 
and a monomial factor, 

p(GA//3)(^) = aou"" exp (-aiu) , (35) 
where m is an integer. The transformation formula ()22j) implies 

These are Gauss-polynomial ensembles including all powers between zero and m. 
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4-5. Gauss-Quartic Ensembles 

We now consider Gaussian probability densities supplemented with a quartic term in 
the exponent, 

P(gQ^)(m) = ao exp (-aiM - Oaw^) . (37) 
With the transformation formula ((221), we obtain 



Q'y^'^pilw) = c-^' exp — —w — —w ^ ■ r , (38) 



2 2 J D_^/2 (ai/v^2, 
where Dp{z) denotes the parabolic cylinder function of order p ^33j. Once more, the 
probability density in superspace contains the functional form of the one in ordinary 
space. However, this example also shows that one can come up with cases in which the 
additional contribution has a rather inconvenient structure. 

4-6. Ensembles Deriving from a Non-extensive Entropy Principle 

An interesting family of ensembles was constructed in Refs. jT^ [T3] . Among other 
features, it yields in a certain parameter range an invariant Levy-type-of ensemble. 
The construction of Refs. [T3] is done for the orthogonal symmetry class, but can 
easily be generalized to all /? = 1, 2, 4. The probability density in ordinary space 

/ K \ 1/(1-9) 

p™(u)=ao(l + ^n) (39) 



A 

depends on a parameter q used in the non-extensive entropy. Moreover, it includes a 
positive parameter k and 

A=^-f, (40) 
q — I 2 

with yU being the number of degrees of freedom (Q). To avoid confusion with the notation 

in the present contribution, we write n, A, /x, instead of a, A, / in Ref. [TH] . 

We consider q > 1. This choice makes the exponent in the probability density ()39p 

negative. Moreover, it requires A > such that 

2 

1 < g < gmax = 1 + - . (41) 
With help of the integral representation ^3] 

oo 

we can obtain the probability density in superspace from the transformation 
formula (j22j) . The u integral has to be done first. Convergence is ensured because 
of the condition A > 0. Collecting everything, we arrive at 



K 



Q™H=c(^)(^l + ^^^j . (43) 

Remarkably, this is again the same functional form as in ordinary space. The (negative) 
exponent 1/(1 — g) in ordinary space is mapped onto —A in superspace. We notice that 
Q^^^^\w) depends on q only via the parameter A which appears twice in Eq. (j43|) . 
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5. Correlation Functions 

After discussing general results for all symmetry classes in Section 15. H we give more 
explicit formulae for the unitary case in Section All results to be given here can also 
be averaged over the external field Hq with some probability density Po{Hq). However, 
as this is an obvious step, we do not go into that further. 

5.1. All Symmetry Classes 

We now have the supersymmetric representation ()14|) for the generating function and 
the one-dimensional transformation formula ()22|1 for the probability density Q^'^^\w) 
in superspace. Hence, we can exploit the advantages of supersymmetry also for norm- 
dependent ensembles. In particular, the level number is, in contrast to the ordinary 
space, an explicit parameter in Eq. ()14|) . Depending on the ensemble, this can makes 
it possible to study the generating function by means of a saddle point approximation 
for large A^. As the details of such a calculation will sensitively depend on the specific 
form of Q^^^^{w), we refrain from attempting a general discussion. 

We can also proceed by observing that Eqs. (fH^jl and pTjl are integral transforms 
of the generating functions involving the Fourier integral. 



+ J) = / f.^ Zf^\x + J, \li1y-) dy , (44) 



J {i2y-Yl^ 

—oo 

or the spread function, 

oo 

{x + J) = J /(^^) (t) (x + J, 2t/(3) dt , (45) 



respectively. Thus, the correlation functions Rl {xi, . . . ,Xk,a, Hq) of all norm- 
dependent ensembles are obtained as single integrals over the corresponding ones 
R^,^^\xi, . . . , Xfc, a, Hq) in the Gaussian case. With Eqs. and (jSj), we find 

i?[.^'^^(xi, . . . ,Xfc,a,ifo) 

— OO 

Rf^\x^pty/(3, Xk^2uJ~^, a, H,^2ty/[3) . (46) 

In those cases where the imaginary unit and the singularities can cause problems, one 
should rather resort to the alternative expression deriving from Eq. ()45|). 

R^^'^\xi, . . .,Xk,a,Ho) 

f(^^\t) 
{2tf/^ 

For Hq = Q and /5 = 1, Eq. (jTTj) was already obtained in Ref. |T7] and a similar result 
is given in Ref. |ll5| for the ensembles deriving from non-extensive entropy. Here, this 
is generalized to all three symmetry classes. 



/ T^^^i^'^'^^^i/^' • • • ' ^'^/^' ^o/V2t) • (47) 
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Of course, supersymmetry is not needed to derive formulae or (|T7j). For Hq = 0, 
one can now use the closed expressions for the correlation functions R^^^^ (xi, . . . , x^, 1, 0) 
from Ref. ^Hj and calculate the correlation functions for all ensembles TOE, TUE 
and TSE. For Hq 0, however, supersymmetry is very helpful, because it provides 
formula ((Zj). As already discussed in Section 1221 the random matrix ensemble, i.e. the 
matrix a, is to some extent decoupled from the external field Hq due to the direct product 
structure. This makes it possible to obtain asymptotic results for large coupling a/D 
as in Refs. |2H1I21I which can then be inserted into formulae and (jUj). 



5.2. Unitary Symmetry Class 

In the unitary symmetry class, i.e. for /5 = 2, we can gain additional insights by following 
the steps outlined in Refs. ^Hl 120] ■ We work with the integral transform 
and the supersymmetric representation of the generating function ((7j). We absorb 
the parameter a into the supermatrix a, which is equivalent to multiplying the 
variance of the Gaussian probability density with a^. The supermatrix a is now 
shifted by X + J to remove these latter matrices from the superdeterminant. We then 
diagonalize a = usu~^ and do the angular integration over the unitary supermatrix u. 
Here, s = diag {su, isu, Ski,isk2) is a diagonal matrix containing the eigenvalues 
■Spi, P = 1, ■ ■ ■ ,k in the bosonic and isp2, p = I, . . . ,k in the fermionic sector. This 
yields [HI 120] 

Zp^(x + J,t) = l-r]{x + J) 

Bk{x + Jj J 

■ exp -trg [s — X — J) 



where the function 

Bk(s) = det 



detg-^[s^®lN-l2k^Ho) , (4^ 



(49) 

p,q=l,...,k 



_Spi - lSq2_ 

is the square root of the Jacobian or Berezinian arising when changing the integration 
variables in superspace to eigenvalue-angle coordinates. The eigenvalues in the 
supermatrix carry a small imaginary increment to ensure convergence, we write s^. 
The function //(x + J) in Eq. ()48p takes care of some Efetov-Wegner-Parisi-Sourlas 
contributions which are not needed here as we are mainly interested in the correlations 
of the imaginary parts. 

Inserting Eq. ()48|) into the integral transform ()45|) and using the normalization (j^UI) . 
we arrive at 

Zp\x + J) = l-r]{x + J) 
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_Le<-. 

detg ^ [s^ ^ In - hk ® Ho) , (50) 
where we introduced the probabihty density 

oo 

e.".W^//'"'Wp^e.p(-|)*. (51) 

Both of the expressions and yield probabihty densities in superspace, the 
former in the fuU, the latter in the eigenvalue space. This is the reason why Eq. ()5ip 
contains the power t'' in the denominator, it arose from the integration over the unitary 
supermatrix. Nevertheless, we can apply the same line of arguing as in Section and 
derive from Eqs. and 1)511) the transformation formula 

Q'^PM = .M/o^ / P^^"\^ + wW^^'-^'^'-'du , (52) 



r((A^2 _2A;)/2) 

with /i = A^^ being the number of degrees of freedom in the unitary case f3 = 2. A 
comparison with the transformation formula ()22|) shows that the probability densities 
in the eigenvalue superspace follows from the one in full superspace by simply replacing 
the number of degrees of freedom A^^ with the reduced number A^^ — 2k where 2k is the 
number of degrees of freedom in the eigenvalue superspace. 



. . 9(Af2-2fc)/2 



(53) 

One also obtains the inversion of Eq. ()52p by modifying the transformation formula (|23|1 
accordingly. 

Hence, we now have an exact expression ()50|) of the generating function for all 
TUE as a 2k dimensional integral with a probability density given by Eqs. ()52|) or (jSHj) . 
Before going over to the correlation functions, a caveat of the same kind as discussed 
in Ref. [201 is in order. The limit a — > can be taken without problems in Eq. (jSUj) . 
thereafter the derivatives with respect to the source variables J can be taken and yield 
the correlation functions for a = 0. Because of some interference with the Efetov- 
Wegner-Parisi-Sourlas term, this should not be done in reversed order for the case 
a = 0. Thus, the following discussion applies to a > 0, where the derivatives of the 
generating function (jMIj) can be taken first. If we may assume that the probability 
density ()52|l does not contain inverses of source variables, we find as in Refs. [TSj [T9\ l20] 

Rp^ (xi, . . . , X,, a. Ho) = ^ / d[s] B,{s) ^Q^P (^trg {s - x)^) 



k N 



«nn^i^. (54) 

p=ln=l ■'pi -'-'On. 

where Hon are the entries of the diagonal matrix Hq. The symbol 53 denotes the 
proper restriction to the imaginary parts which will be explained below. As already 
observed in Ref. ^H]; the determinantal structure of the GUE correlation functions 



Norm- dependent Random Matrix Ensembles in External Field and Supersymmetry 15 



arises in this supersymmetry approach as an immediate consequence of the determinant 
structure (jlH))- Thus, it follows form the Jacobian in superspace. The only term in the 
integrand which can destroy this feature for the TUE is the probability density (w). 
To circumvent this problem, we use formula ()51|) which makes possible to advantage of 
explicit results for the GUE correlation functions in the presence of an external field. 
The formulae given in Ref. 120] are slightly more general than what we need in 
the present context, because they also contains an integral over the probability density 
Po{Ho) of the external field Hq. However, a 6 function Po{Ho) trivially yields 

Rf^\xi,. . . ,Xk,ta'^,Ho) = det\CN{xp,Xg,ta'^,Ho)] . (55) 

The kernel is given as the double integral 

CN{xp,Xg,ta'^,Ho) 

1 dsids2 ({is2 — Xqf' isi — Xpf\ 

exp 



27r2ta2 J J si-is2 \ 2ta^ 2ta^ 



-oo —oo 



3 



n "^^^ , (56) 



^ 182 - Ho 
n=l ^1 ~ Hon 

where we drop the indices p and q of the integration variables. The present notation 
slightly deviates, in a hopefully self-explanatory way, from the previous one. We now 
have the variance ta^ as an argument, because it contains the parameter a after the 
above mentioned changes of integration variables. The process of going over to the 
imaginary parts of the correlation functions amounts to inserting the definition 

TT ^■^S — 1 f JL iS2 — Hon iS2 — -ffon \ ,r-^\ 

^[[^ 73- = ^ M1t= 73--11t+ 73- ■ (57) 

n=l ■'l ~ -"On \„-i Si — SlQn ^=1 ^1 ~ -"On / 

We notice that the term l/(si — is2) in the integrand of Eq. (j56|) is the remainder of the 
Jacobian. After Refs. [IHl UHl 120] , such double integral expressions were also derived by 
other authors. 

We combine these findings and arrive at 

R^P\xi, . . .,Xk,a,Ho) 

= dt——- det CN{xp,Xq,ta'^, Ho) (58) 

^ ^ 

which is an exact representation of the TUE correlation functions for finite level number 
in the presence of an external field. For convenience, we derived this result using the 
spread function f^'^'^\t). However, depending on the specific form of the probability 
density Q^^'^\w) given in Eqs. ()52j) and (jK!^ . one might want to prefer another integral 
representation. Any such representation which involves a Gaussian will lead to a result 
of the form ()58|1 . but with the spread function replaced by another function. It should 
also be mentioned that the result ()58p can alternatively be derived starting directly 
from Eq. (jlTj). This, however, would not lead to the probability density Q^^'^\w) and 
the corresponding observation concerning the reduced number of degrees of freedom. 
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Furthermore, we notice that the double integral (|56p for the kernel can be evaluated in 
closed form. For the sake of completeness, we give this result in Appendix B 



6. Summary and Conclusions 

We showed that the all norm-dependent random matrix ensembles TOE, TUE and 
TSE have a supersymmetric representation. Hence, supersymmetry is by no means 
restricted to Gaussian probability densities. We mapped the functions generating the k 
point correlation functions onto their supersymmetric analogues. All our results include 
an external field. No approximation was made, all expressions are exact and for finite 
level number. We derived transformation formulae which yield the probability density 
in superspace as one-dimensional integral involving the probability density in ordinary 
space. These formulae state the most important conceptual insight of this contribution. 
We emphasize that only the existence, but not the explicit knowledge of the Fourier 
integral or the spread function is needed to obtain the probabihty density in superspace. 
The transformation formulae clarify the mechanism of how supersymmetry manages to 
reduce the number of degrees of freedom. We worked out several examples. Remarkably, 
the functional forms of the probability densities tends to be very similar in ordinary and 
superspace. In particular, this is so for the whole family of ensembles deriving from a 
non-extensive entropy principle. 

From a practical viewpoint, our most important findings are a series of exact results 
for the correlation functions which can be used in applications. Different limits for the 
level number or other parameters can now be studied depending on the ensemble of 
interest . The most explicit formulae are for the TUE where we employed the determinant 
structure of the GUE correlations. We also derived a probability density for the TUE 
in the superspace of eigenvalues. 

Can supersymmetry be applied to ensembles which are more general than the 
norm-dependent ones? — Indeed, a supersymmetric representation is possible under 
quite general conditions. However, as this construction requires a completely different 
approach, we defer it to another contribution 34j. 
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Appendix A. Moments of Norm— dependent Probability Densities 

The following calculation, although rather straightforward, yields an interesting side 
result which might be useful in other applications. This is why we sketch the calculation 
here. Inserting the diagonalizations H = UxU^^ with x = diag (xi, . . . , xtv) and a 
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double degeneracy for /3 = 4, we find from the definition of the moments in Eqs. (jl2|) 

j\^(T/3) ^ ^m{N-i)/4_ r^(i + (3/2) 



nlir(l + n/3/2) 
J p(^^)(Trx2) (Trx^ |Ajv(x)|^rf[x] , (A.l) 

where Aj^{x) = Y[n<m{^n — a^m) is the Vandermonde determinant. The constant in 
front of the eigenvalue integral contains the result of the integration over U and also 
some factors stemming from the Jacobian of the transformation to eigenvalue-angle 
coordinates. We view the eigenvalues as components of a vector r = (xi, . . . ,xn) in N 
dimensions and introduce polar coordinates f = re where r is the length and e a unit 
vector depending on — 1 angles. The volume element reads d[x] = d^x = r^~^drdVL 
where dVL is the infinitesimal sohd angle. Because of = Trx^, we have 



J^(TI3) ^ ^f3N(N-l)/i 



r^(i + /?/2) 

n;:=ir(i + n/5/2) 

oo 

A^(e)|^dl] /'r'^+^-i+/3^(^-i)/2p(^/3)(r2)cir . (A.2) 





The angular integral can be infered by inserting the Gaussian probabihty density (0)) 
and considering i/ = 0, i.e. Mq"^'' = 1. The radial integral can then be done and we find 
as an interesting side result 

\^^^^)\ - 2miv-i)/4-ir^(i + ^/2)r(^/2) • ^^-^^ 

Putting u = r"^ and collecting everything we arrive at the second of Eqs. p2j) . 
Appendix B. Evaluation of the Kernel 

We start by observing that the determinant can be written in the form 

N ■ IT N / ■ \ 



-1 S'l ~ \ Si — Hon / 

„=1 Si - Hon m=Q ll^,gf^W l^On " ^Om'j 

Here, f^^^ is a subset of the A^ — 1 indices remaining when removing the index n from 
the original A^ indices such that m pairs (n, m') are formed. This can also be formulated 
in terms of symmetric functions. For example, in the case A^ = 3 and tt, = 1, one has 

^ {lS2 - Si)" 



Om'J 

2 



^ ^ ^S2 - Si ^ iS2 - Si ^ [iS2 - Si) 2^ 

Hqi — Hq2 Hqi — Hq3 {Hqi — Hq2){Hqi — Hq^) 
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This yields a most convenient expression containing only 6 functions, 

^ is2 - H ^ 
S n ~ 77^ = vr ^(is2 - si)6{si - Hon) 

n=l ~ -"On n=l 

which facilitates the evaluation of the si integral in Eq. (|56|l . Importantly, the difference 
Si — is2 also disappears in the denominator and the remaining S2 integration simply 
becomes 

\2' 



J {iS2 - Hon) exp I I dS2 



with lirn{z) denoting the Hermite polynomial of degree m. Collecting everything, we 
find 

^ , , 2 rr, ftC^'Y^" ^ exp((iJo„-X,)V2to^) 

Cn{xp, Xq,ta ,Ho) = \ — \ }^ ^=== 

H„((Xg-i7on)/V2to2) 



We notice that the first sum extends from zero to — 1, exactly as in the case of the 
GUE without external field. 
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